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ABSTRACT 

The results of two different methods for determining the 
rigidity modulus of a soft sediment are compared. In one method 
the resonant characteristics of a torsionally oscillating rod 
which are sensitive to the shear acoustic impedance of the sedi- 
ment in which the rod is imbedded determine the complex rigidity. 

The second method utilizes the observation of the phase velocity 
of an interface wave at the water-sediment boundry. Shear wave 
speeds computed from the experimental data from both methods are 
quite similar in magnitude. For the sediment used here, the 
average value of shear wave speed determined from the interface 
wave experi.ment was 33 m/sec while the shear wave speed determined 
from the measured rigidity was 29 m/sec. The difference lies with- 
in experimental uncertainty. Trends in the mass-physical properl ies 
of the sediments are investigated by comparing graphically the de- 
pendence of both the real and imaginary parts of the complex rigid- 
ity on density, porosity, sound speed, silt and clay percentages, 
Poisson ^s ratio and density times sound speed squared. 
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I. 



INTRODUCTION 



Growing concern for the accurate and comprehensive description 
of the acoustic reflection and attenuation characteristics of the 
ocean bottom is continually stimulated by attempts to fill the 
large gap in predict ibility of the actual reflection losses* The 
development of a model for the reflection coefficient which more 
accurately describes the physics of transmission and reflection 
should result in improved capability for predicting transmission 
loss for sound in the ocean in those cases where bottom interaction 
is important, such as passive and active sonars in shallow water 
and active sonars using the bottom bounce mode. 

Historically, the ocean bottom has most often been modeled as 
an interface between two fluids of differing densities and sound 
speeds. Usually referred to as the Rayleigh model, it is generally 
accepted because normal* surficial sediments have high porosities 
and low rigidities which approximate those of a fluid. The 
Rayleigh model is sometimes found to be deficient, especially at 
lower grazing angles [ 3 , 9 ]- 

Some have accounted for this disparity by presuming that the 
sediment posseses a finite rigidity and hence is capable of sup- 
porting -a shear wave [3]* Using a viscoelastic model and complex 
Lame* constants, one is capable of describing the formation of 
shear-interface waves and absorption of compressive sound waves 
in the sediment. Bucker , et al. [3] finds better agreement with 
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experimental measurements for the theoretically calculated bottom 
loss using this viscoelastic model and associated sediment para- 
meters obtained by Hamilton [9jl0,ll] than for the Rayleigh model 
in a number of cases. 

In order to obtain both complex Lame* constants for the vis- 
coelastic model it is required that, in addition to the measure- 
ment of sound absorption, both shear wave speed and compress ional 
wave speed be measured. Hutchins [ 12 ] demonstrated the feasibility 
of measuring the complex shear modulus of soft sediments (artifi- 
cial) by noting the effects caused by the rigidity on the resonant 
behavior of a torsionally oscillating rod imbedded in the sediment. 
Utilizing a laboratory prepared kaolin-water mixture, he developed 
the viscoelastometer technique in the frequency range of 3 S .3 io 
38.9 KHz. 

Cohen [5lj by improving the design of the instrument, increased 
the range of operation to 5.8 to 38 KHz. Cohen used both kaolin- 
water and bentonite-water combinations as artificial sediments. 

Utilizing Hutchins* and Cohen *s observations that soft, 
clayey, simulated sediments did, in fact, have a measureable shear 
modulus, Bieda [ 2 ] continued this research by developing a new 
viscoelastometer suitable for use in making in situ measurements 
and the subsequent laboratory analysis of 20 sediment cores ob- 
tained in the Monterey Bay area. 

The results obtained from use of these probes [2,5,12] are 
quite encouraging and are certainly within reasonable bounds for 
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this type of sedimentary material (E.L. Hamilton, personal 
communication). However, there has been no independent deter- 
mination of the rigidity for comparison and verification of the 
method. Consequently it was required that experiments be con- 
ducted that would test the validity of the probe method. Suc- 
cessful development of such a probe would make the necessary 
measurements of ocean bottom parameters a more convenient 
operation . 

Another method for determining the elastic constants involves 
measuring directly both the compressional wave and shear wave 
speeds; knowing these, it is a simple matter to determine the real 
elastic constants. The imaginary parts of the elastic moduli must 
be obtained from measurements of the attenuation of the compress- 
ional and shear waves. The former method is used in this experi- 
ment as an independent means for validating the measurement of 
the real components of the rigidity made by the viscoelas tometer 
developed by Bieda [ 2 ]. The method described is essentially the 
same as that used by Hamilton, et al. [ll] in measurement of the 
interface wave speed and compressional wave attenuation in sedi- 
ments in the San Diego area^. 



Some controversy appears to exist as to just what title should 
be applied to the interface wave that exists at the suface between 
an elastic, solid bottom and the overlying water. Although gener- 
ally referred to as Stoneley waves, it appears that Stoneley [l?] 
never included the possibility that the upper medium could be a 
fluid. Other people tend to choose the name Scholte waves since 
Scholte [l6j generalized the nomenclature to include fluid-solid 
interface waves. Since the environment in which the author ^s data 
were taken did not entirely approach an elastic-solid bottom of 
infinite extent below a infinite half-space of homogeneous fluid, 
it was decided to use the term “interface wave“ when describing the 
elastic wave existing at the interface between the bottom and the 
overlying water. 
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II . EQUIPMEm: AND PROCEDURES 



A. LOCATION OF THE EXPERIMENT 

Cohen [ 5 ] has shown the viscoelast ometer to be ineffective in 
sediments of large mean grain size, possibly because the grain 
dimension becomes comparable to the shear wavelength in the sandy 
sediments and therefore the assumption of homogeneity becomes 
invalid. Thus there was a need for finding a sediment which was 
silty clay or clayey silt. 

There are vast areas of mud bottoms in the Monterey Submarine 
Canyon area of Monterey Bay, California, that would normally 
suffice for this type of experiment. However, since the services 
of a submersible platform were not available the evaluation had 
to be made at a quite shallow depths (1 to 1.5 m) near shore. A 
suitable mud sediment site was found in the tidal flats of the 
Kirby Park area of Elkhorn Slough, a few miles north of the city 
of Monterey. 

B. INTERFACE WAVE MEASUREMENT TECHNIQUE 

Direct measurement of the speed of shear waves in very soft 
marine sediments is made virtually impossible at acoustic fre- 
quencies because of the rapid attenuation of the shear waves and 
the difficulty of coupling a suitable transducer to the medium. 

It has, however, been determined by Strict and Ginzbarg [18] and 
by Hamilton, et al. [ll] that the ratio of the shear wave speed 
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to the interface wave speed remains relatively constant for a 
given value of the ratio of the two densities of adjoining media 
irrespective of other varying bottom constants. 

While compressional wave speed and wet density of the sediment 
can easily be measured in the laboratory and corrected to in situ 
values, it was required that the interface wave speed be measured 
and corrected to shear speed as previously accomplished by 
Hamilton, et al. Cll]. 

An interface wave, generated by the detonation of a blasting 
cap at the water -sediment interface, is examined using three or 
four vertically-polarized, velocity-sensitive geophones mounted at 
known intervals on the interface. An aluminum rod and flexible 
lines served as a convenient support and the means for spacing of 
the geophones and the electrically-fired blasting caps which were 
mounted under a heavy brass plate. 

During the firing of these blasting caps, the response of the 
geophones is monitered at a remote station. The signal from the 
geophones passes through an impedance-matching transformer, a set 
of attenuator pads and a low-noise amplifier capable of as much as 
50 db gain to a Honeywell Visicorder where the subsequent amplitude 
versus time responses are plotted as shown in Figure 2. A sketch 
of the operational assembly is shown in Figure 1. 

The Vis icor der traces are studied to determine the points of 
constant phase in each of the geophone traces (point a in Figure 2) . 
This is generally best accomplished by attempting to match waveforms 
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which is somewhat aided by irregular spacing of the geophones and 
adjustment of gain so that amplitudes are comparable. Having 
found the points of constant phase, and using convenient time 
reference marks, the interface wave speed is easily determined. 

Having obtained the interface wave speed, the graphs provided 
by Strict and Ginzbarg [l8] are entered with the ratio of the 
density of the overlying sea-water to the density of the sediment. 
The graphs are plotted with the ratio of interface wave speed to 
shear wave speed versus the ratio of shear wave speed to compress- 
ional wave speed, varying the par ameter rat io of water density to 
sediment density. Since both the ordinate and abscissa are ratios 
involving the unknown, shear wave speed, an additional assumption 
about the magnitude of shear wave speed allows one to accurately 
determine the necessary value. As the ratio of shear wave speed 
to compr essional wave speed approaches 0.5, the individual graphs 
all begin to approach a horizontal asymptote which may conveniently 
be read from the ordinate with a maximum possible error of 
5 percent . 

Having determined the value of the ratio of the interface wave 
speed to shear wave speed, and knov/ing the interface wave speed, 
it is easy to compute the shear wave speed (V^ = V^/A, where A is 
the ratio, the shear wave speed and the interface wave 
speed) . ' 
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C. THE VISCOELASTOMETER 



1 • Equ ipment 

The apparatus designed by Bieda [ 2 ] is a result of a 
series of improvements on the original devices constructed and 
evaluated by Cohen [ 5 ] and Hutchins [ 12 ]. It consists of a 
10^- inch long, %-inch diameter cylindrical rod which is excited 
in torsional oscillation by a centrally mounted barium titanate 
transducer. The rods, constructed from a constant modulus alloy, 
Nl-Span-C, were heatt^eated and highly polished. The electrical 
leads to the transducer are passed through a supporting yoke which 
supports a clamp that holds the transducer in position. The clamp 
is secured to the transducer at its center, a node of torsional 
motion for all resonant modes in the rod. A sketch of the 
apparatus is shown in Figure 3* Various compounds were used to 
maintain the water-tight integrity of the electrical connections, 
generally RTV insulating rubber was deposited in the large crevices 
at the junction of the clamp and the transducer while the other 
boundries were coated with varnish or glyptol. 

A .Dranetz lOOB Impedance- Admittance meter , a suitable 
oscillator and a frequency counter were used to measure the elec- 
trical conductance and the resonant frequencies. 

2 . Theory and Method of Measurement 

When the viscoelas tometer is excited at its mechanical 
resonance, a torsionally oscillating standing wave exists along 
the axis of the rods and transducer. The torsional motion of the 
walls of the rod generates a shear wave which propagates radially 
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outward into the sediment. The radiation impedance of the sur- 
rounding sediment has a modifying effect on the magnitude of the 
electrical resistance of the probe and its frequency of mechanical 
resonance, compared to those values measured in air (approximately 
no-load conditions). The change in resistance and resonant fre- 
quency are thus direct functions of the rigidity of the sediment [l4]. 

Under the assumption that the shear waves that propagate 
outward from the viscoelastometer have wavelengths much smaller than 
the radius of curvature of the rods and that the amplitude of these 
shear waves is very rapidly attenuated, the waves propagated in the 
sediment may be considered to be plane waves with the formula: 

T = exp [ \j rr:^/r] (l+i)z 3 (1) 

where T is the shear stress amplitude, T^ the initial amplitude, 
f the frequency, and r? the shear viscosity. The harmonic time- 
varying function has been suppressed for sake of convenience. 

In a Newtonian fluid, the coefficient of shear viscosity 
• . 
is ?7 = T/S, where S is the shear strain and S = bs/5t. In our 

viscosity model a complex rigidity (G^) l^^is been assumed along 

with a complex viscosity {V * Thus, 

T?c = ^1 - 

In a Newtonian fluid, 77 ^ = 0, and = 0 in a perfectly elastic 
solid. Since 77 ^ = T/iCOS = - iG^/ CO for simple harmonic motion, 
where CO is the angular velocity, 

and . (3) 
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The specific radiation impedance (Z) of the sediment is 



the ratio of shear stress to particle velocity. Following 
Mason [l4l and McSkimin [iSl? for a Newtonian fluid 

Z = R + iX = TTf ?7p (1 + i) (4) 



where R and X are the specific acoustic impedance and reactance 
of the sediment, respectively, and p is the sediment density. For 
the case when viscosity is complex, the following relations result 

= ^ T7 = 

^1 COO ^2 COO 

(5) 



G 



1 



2 2 
R -X 



2RX 
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Thus, it becomes apparent that determination of the complex 
rigidity or viscosity requires only a knowledge of real and 
imaginary parts of the specific acoustic impedance for plane 
shear waves. 

Mason [l4] has shown that constants of proportionality can 

2 

be determined which relate the total acoustic resistance (R^) to 
the change of electrical resistance ( AR ), i.e., AR = K^R^, and 
the total acoustic reactance (X^) to the change of resonant fre- 
quency (Af), i.e., Af = K^X^ , when the loading on the rod is 
changed from air to the medium under study. Here, R = ^^id 

X = X^/1, where 1 is the length of the probe immersed in the 
sediment. Thus, measurement of the acoustic impedance is reduced 
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Transducer 



which drives the rod. 
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to measurement of the change observed in the frequency of 
resonance and in the resonant value of electrical conductance 
when the probe is operated in the material to be analyzed. 

Evaluation of the constants and is accomplished by 

measurement of AR and Af in Newtonian calibration fluids. Since 

e 

T}^ = 0 in a Newtonian fluid, R = X; knowing the viscosity and 
density of the test fluid, the procedure for calculating and 
is straightforward. 

3 • Core Analysis 

In the same relative location of the interface wave experi- 
ment ten cores were taken for purposes of v iscoelastometer analysis, 
density and compress ional wave speed measurement. While maintaining 
the core at its original temperature, the complex rigidity is de- 
termined using Bieda’s equipment. The measured real part of rigid- 
ity and the observed density are used to determine a shear wave 
speed, i.e., = (G^/p) . 

The core samples were analyzed at two of the resonant 
frequencies, 27 kHz and 38 kHz, due to anomalous behavior at the 
lower frequencies, especially in the stiffer samples. Since the 
constants and tend to vary with frequency and viscosity, an 
average value was used at each frequency from calibrations in 
fluids with viscosities roughly equivalent to those expected in 
the sediments. Table III contains examples of the variance of 
these constants over temperature and viscosity. 
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The computed shear velocities are averaged and are com- 
pared to the value for shear velocity obtained in the interface 
wave experiment in Table II. 

D. COMPRESSIONAL WAVE SPEED MEASUREMEI>IT 

Compress ional wave speed is measured in each of the sample 
cores by determination of the time required for a 2 MHz com- 
pressional pulse to pass through a known length of sediment. 

Barium titanate disk transducers are used as a transmitter and 
as a receiver , while the transmitted pulse width and frequency 
are varied in close proximity to 2 MHz in order to obtain resonant 
operation of the transducers. 

The two transducers are mounted facing each other as shown 
in Figure k* The sample, approximately 3 inches long, is inserted 
in the holder which is in contact with one transducer; the other 
transducer can then be moved to contact the other end of the 
sample. A dial indicator is used to measure the separation while 
a time mark generator and the output of a triggered pulsed oscillator 
are displayed on the face of a dual trace oscilloscope. The tem- 
perature is maintained at a constant value, identical to that 
measured during the interface wave experiment. 

Calibration of the velocimeter is accomplished by measuring 
a sequence of time delays for various transducer spacings in a 
mixture of 22 parts ethanol to 100 parts distilled water. This 
combination was shown to be relatively insensitive to small tem- 
perature changes in the neighborhood of the temperatures of this 
experiment [2l]. 
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The sediment sound speeds and the speed of sound in the water 
samples (computed using Wilson’s equation) taken in the vicinity 
of the experiments are shown in Table I. 

E. WET DENSITY AND POROSITY 

Wet density is measured by weighing a known volume of the 
saturated sediment. The weight of this same sample after heating 
in an oven at 105^0 for a period of 24 hours is again determined. 
If it is assumed that the weight lost is completely derived from 
evaporation of the entrapped water and that the weight of 1 mX of 
water is Ig , the porosity (P) is easily calculated, 

_ Volume of Voids 

Total Volume ’ ^ 



F, CALCULATED PARAMETERS 



The shear wave speed is calculated by use of the real part of 
the measured rigidity, , and the sediment density, p, in the 
for mula : 




( 7 ) 



Shear wave speed is also computed from inteface wave speed by the 
formula : 



Vg = V^/A , (8) 

where A is the constant of proportionality derived from the graphs 
of Strict and Ginzbarg Cl8j. Poisson’s ratio Cr, is calculated by 
the relation: r 1 



-I 



((Vp/V^) -1) 



( 9 ) 



19 



k 




where is the compress ional wave speed and is the shear wave 
speed. The measured and computed physical constants are shown 
in Table I . 

G. SAND- SILT -CLAY RATIOS: MINERALOGY 

Sand-silt-clay ratios were calculated at the Naval Postgraduate 
School by the sieve pipette method. Clay mineralogy analyses, 
performed by Dr. W.R. Bryant, Texas A&M University, are shown in 
the Appendix. 
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III. LIMITATION OF THE METHODS 



A. RESONANT VISCOELASTOMETER TECHNIQUE 

While usage of this probe is not strictly limited to the 
laboratory, a very strong dependence on temperature of the con- 
stants and (Table III) requires that measurements for com- 
parison be made in a temperature-controlled environment. To 
allow for correct measurements to be made from 2-3/4 inch diameter 
cores, the assumption is made that the shear wave length is much 
smaller than the distance from the probe to the core liner. It 
was noted that at the lower frequency, especially in stiff er 
sediments, erratic readings were generally obtained, possibly from 
the violation of this assumption. This limitation should not pose 
any great problems when the probe is used in situ as there should 
then be no surrounding barrier in coupling to the motion of the 
oscillating rods. 

It was also noted that the values of and vary with 
changing viscosity as well as frequency and temperature. Con- 
sequently, the final constants used were obtained from an average 
of those values obtained over the lower values of viscosity 
(1-15 poise). This could cause an error in of no more than 
30 percent and perhaps 110 percent in the value of • 

Another obvious limitation to the laboratory analysis is that 
no matter how carefully a core is taken, there is always some 
structural distortion of the sample as well as some drying before 
analysis can be made. 
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B. INTERFACE WAVE TECHNIQUE 



The most serious breach of assumptions occurs in the measure- 
ment of the interface waves at shallow depths. The determination 
of the interface wave speeds are all predicted on an interface 
between liquid and elastic-solid half spaces [4,7l. In the mea- 
surement of the lower frequency (about 10 Hz) interface waves in 
shallow water it would seem that the effect of the overlying water 
is quite small, since the wavelength for 10 Hz compr essional waves 
in water is so large compared to the depth. Thus, the situation 
is very nearly that normally encountered for Rayleigh waves. This 
effect could cause the ratio of the interface wave speed to the 
shear wave speed, A, to approach a value of 0.92 C4l. This could 
yield a possible error in the measured shear wave speed of as much 
as 10 percent . 

Higher frequency compressional waves, although much faster, 
occasionally tend to mask the arrival of the interface wave due to 
the exponential tail of the explosion. Often this masking caused 
some uncertainty as to the actual arrival time of the interface 
wave.. The standard deviation of the measured values of interface 
wave speed are about 11 percent of the average measured value. 

Since the distance between the geophones enters directly into 
the computation of the interface wave speed, the uncertainty of 
about 2 'em could yield an error of about 5 percent in the result. 

Since only quantitative information describing the relation- 
ship between and interface wave an a solid shear wave is that 
given by Strict and Ginzbarg [l8], the author was forced to make 
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the approximation that shear wave speed would always be less than 
one-half the compr essional wave speed in order to obtain a value 
for the needed ratio, A, as noted in section II-B. This uncer- 
tainty in the constant, A, could yield an error of no greater 
magnitude than 5 percent. 
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IV. RESULTS AND DISCUSSION 



Table II presents a comparison of the two methods of measuring 
rigidity by showing the average values of the shear wave speed as 
calculated from measurements of the real part of the rigidity and 
the shear wave speed computed from the interface wave experiment. 
Extremely good correlation between the two values tencfe to confirm 
the validity of the torsional viscoelastometer results. tt 
must be pointed out, however, that the viscoelastometer was 
operated at frequencies of 27 kHz and 3 S kHz while the interface 
wave frequencies were about 10 Hz. Although there is a great 
disparity in frequencies, there is no reason to expect that 
dispersion of shear speed should occur [ll]. The viscoelastometer 
technique is subject to a maximum precision of 30 percent in 
calculating the real part of the rigidity, largely arising from 
an inability to accurately measure all the parameters. While 
these uncertainties of the method point toward the need for 
stricter maintenance of conditions during viscoelastometer mea- 
surements and during calibration, they are- not of sufficient 
magnitude to discredit the method. 

The mechanical disturbance in coring and in transporting the 
samples. very likely disturbs the physical properties by an inde- 
terminant amount. It is also known from a visual inspection of 
the cores that the bottom in this area was far from being homo- 
geneous, in either the horizontal or the vertical plane. 
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The measured physical parameters for the ten cores taken tend 
to follow the same general patterns (Figures 5 to 19) as did those 
presented by Bieda [ 2 ] with the major exception that the cores 
taken in the slough possess less than 7 percent sand, making 
them much more acceptable for viscoelastometer analysis. 

It is significant to note that the rigidities obtained ranged 

6 62 

between 5.1 x 10 to 2.7 x 10 dynes/cm for tbe real part and 

6 62 

1.1 X 10 to 7.9 X 10 dynes/cm for the imaginary part. These 
values, as pointed out by Bieda [ 2 ], are significantly lower than 
those found by Hamilton [ll], but additional evidence continues 
to point toward the direction of lower rigidities in surficial 
sediments [l,6]. 

The trends noted by Bieda were also observed in this research; 

the real part of rigidity increases with increasing density, 

2 

increasing percent silt and with increasing pv^ . Rigidity de- 
creases with increasing porosity and percent clay. The same 
trends are noticeable but less marked with the imaginary part due 
to the greater degree of scatter caused by experimental imprecision. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



The relatively close agreement of the two values of shear 
wave speed strongly tends to validate the usefulness of the 
viscoelast ometer technique for measuring the real part of the 
complex rigidity in soft, fine grained surficial marine sediments. 

The mass-phys ical trends noted by Bieda [2] are seemingly 
substantiated by these observations but require a great deal more 
study to completely classify definite behavioral functions. 

In order to further exploit the use of the viscoelastometer , 
it is recommended that further efforts be applied to the design 
of the viscoelastometer in order to minimize the effects on its 
calibration of grain size, stiffness and temperature. 

Further effort is also indicated toward improving the precision 
in measuring the imaginary part of the complex rigidity and toward 
its validation, utilizing a measurement of the sound absorption in 
the sediment as a test. 

Construction of an artificially prepared ocean bottom (kaolin- 
water) is recommended in order to circumvent the damage done to the 
cores when taken from the bottom and to improve homogeneity and 
temperature stability of the medium. 
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TABLE I. Tabulated Physical Properties 



SAMPLE 


SOUND 

SPEED 

m/sGC 


SHEAR 

SPEED 

m/sec 


1. 


1487.0 


22.1 


2. 


1502.0 


21.5 


3 . 


****** 


32.9 


4 . 


1510.0 


21.1 


5 . 


1505.5 


23.5 


6. 


1506.1 


26.4 


7 . 


1517.3 


47.0 


8. 


1510.1 


33.5 


9 . 


1493.8 


33.1 


10. 


1520.1 


*28.8 



CSED/ 

CWAT POISSON'S RATIO 



0.983 


0.4998, 


0.993 




***** 


******* 


0.998 


0-'*99,o 


0.995 


0.49987 


0.995 


«-W9gi, 


1.000 


0.499.J 


0.998 


0.499yj 


0.987 


0.49975 


1 .000 


0.49982 



****** Sound speed unattainable 

* Compt essional speed in sediment/compr essional speed in ivater 
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TABLE I. Tabulated Physical Properties (cont.) 



SAMPLE 


G -1 

dynes/cm 


G -2 

2 

dynes/cm 


DENSITY 

g/cm^ 


Eorosity 

% 


1 . 


5.8 


2.5 


1.21 


65.9 


2 . 


5.1 


3.0 


1.25 


64.7 


3. 


14.8 


6.9 


1.35 


74.4 


4 . 


6 .0 


1.8 


1.36 


78.3 


5 . 


7.6 


1.1 


1.37 


72.9 


6 : 


9.8 


2.7 


1.40 


74.4 


7 . 


27.8 


7.9 


1.25 


72.5 


8 . 


17.3 


7.0 


1.54 


68.9 


9 . 


16.1 


5.2 


1.46 


71.1 


10 . 


12.2 


5.8 


1.47 


72.0 
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TABLE I. Tabulated Physical Properties (cont.) 



AMPLE 


PERCENT 


PERCENT 


PERCENT 


UMBER 


SAND 


SILT 


CLAY 


1. 


2 . 8 % 


29.4% 


67.8% 


2. 


7.2% 


52.6% 


40.4% 


3. 


2.2% 


70.9% 


26.9% 


4. 


6.5% 


25.7% 


67.8% 


5. 


1.4% 


56.5% 


42.1% 


6. 


1.9% 


70.8% 


27.3% 


7. 


2.0% 


72 . 2 % 


25.8% 


8. 


1.9% 


79.6% 


l8.5% 


9. 


2.2% 


79.7% 


l8.l% 


10. 


2.2% 


72.3% 


25.5% 
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TABLE I. Tabulated Physical Properties (cont.) 



**** 



SAMPLE QC? 

2 , 3 2 . 

(10 g-m /cm -sec ) 



1. 


2.67 


2. 


2.82 


3. 


**** 


4. 


3.67 


5. 


3.10 


6. 


3.17 


7. 


2.87 


8. 


3.52 


9. 


3.26 


10. 


3.40 



Sound speed unattainable 
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TABLE II. Comparison of Shear Speed Results 



AVERAGE STANDARD DEVIATION ESTIMATE OF 

VALUE OF DATA POINTS PRECISION 



Ratio of V./V 
1 s 


0.758 


0.062 


Shear wave speed 
(from measured 
r igidit ies ) 


29.0 m/sec 


8 m/sec 


Interface v;ave speed 
(from interface wave 
exper iment ) 


28.8 m/sec 


4 m/sec 


Shear wave speed 
(from above) 


33-3 m/sec 


4 ra/sec 



+ 5 % 

+ 30% 



+ 10 % 



+ 15% 
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TABLE III 



Variation of Constants After Averaging at Low Viscosities 



Temperature ( °C) 


4 kHz 


27 kHz 


38 kHz 


10 


6.04 


0.43 


0.34 


17 


5.73 


0.69 


0.62 


22 


6.20 


0.71 


0.70 






X 10“^) 




Temper a tur e ( 


4 kHz 


27 kHz 


38 kHz 


10 ■ 


1.36 


1.38 


1.52 


17 


2.96 


2.67 


2.83 


22 


3.81 


2.66 


2.72 
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Figure 1. 



Sketch 



of Interface Wave Measurement 



Apparatus 
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d=1.5 m 




Figure 2. Visicorder Prints from Interface Wave 
Measur ement 
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INTERFACE WAVES 
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Figure 3* Sketch of Viscoelast oraeter and Associated 
Equipment [ 2 ] 
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Figure 4. Sketch of Compress ional Wave Velocimeter 
Apparatus [ 2 ] 
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POROSITY (%) 

Figure 6. as a Function of Porosity 
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Figure 7- as a Function of Sound Speed 
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POISSON'S RATIO 

Figure 8. as a Function of Poisson’s Ratio 
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Figure 9. as a Function of Density and Sound Speed 
Squared 
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Figure 10. as a Function of Percent 


Silt 
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Figure 11. as a Function of Percent Clay 
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Figure 12. as a Function of Density 
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Figure 13* a Function of Porosity 
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Figure 15. as a Function of Poisson's Ratio 
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Figure 19 • Shear Speed as a Function of Sound Speed 
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APPENDIX: CLAY MINERALOGY 



The results of the quantitative analyses of the clay minerals 
from X-ray diffraction are shown in Table IV. The relatively high 
percentage of chlorite in these samples is due to river drainage 
from areas of outcrops of the Franciscan Formation. The graphs of 
G and G^ as a function of the various clay minerals show the data 
from this report as well as those from Bieda ( 1970) ( Figures 20 - 25). 
Little in the way of definite trends may be seen, possibly due to 
the fact that these data are all from nearshore sediment samples. 
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TABLE IV. 



Clay Mineralogy 





PERCENT 


PERCENT 


PERCENT 


PERCENT 


SAMPLE 


MONTMORILLONITE 


ILLITE 


KAOLINITE 


CHLORITE 


2. 


72.88 


13.56 


6.51 


6.98 


3 . 


69.31 


13.97 


8.05 


8.66 


4 . 


60.72 


17.86 


11.82 


9.60 


5 . 


51.35 


21.62 


12.16 


14.86 


6. 


47.33 


25.95 


11.76 


14.96 


8 . 


64.91 


12.21 


10.53 


12.37 


9 . 


48.82 


24.24 


12.53 


l 4 . 4 l 
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% Montmorlllonite 

Figure 20. as a Function of Percent Montraor illonite 
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Figure 21. as a Function of Percent Illite 
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Figure 22. as a Function of Percent Chlorite 
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Figure 25- G 2 as a Function of Percent Chlorite 
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